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Abstract. The intra-day variable source 0917+624 dis- 
plays annual changes in its timescale of variability. This 
is explained in terms of a scintillation model in which 
changes in the variability timescale are due to changes in 
the relative velocity of scintillation pattern as the Earth 
orbits the sun. 
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1. Introduction 

The radio source S5 0917+624, identified with a z = 1.446 
quasar (Stickel & Kuhr 1993), was one of the first discov- 
ered to show strong intra-day variability (IDV) at radio 
wavelengths (Heeschen et al. 1987). Recently, Kraus et 
al. (1999) reported three-epoch 5 GHz observations, De- 
cember 1997, September 1998 and February 1999, which 
show significant changes in the IDV characteristics of 
this source. They suggest several possible interpretations 
which have at their basis structural changes, and hence 
brightness temperature changes, within the source itself. 

Alternatively we suggest that the observed changes 
in IDV result from variations in the interstellar medium 
(ISM) velocity as the Earth revolves around the sun. We 
present evidence to support this suggestion based on pub- 
Hshed data on the IDV of 0917+624. 



2. Changes in the Intra-Day Variability in 
0917+624 

Kraus et al. (1999) undertook to explain the observed 
changes in the IDV properties as an intrinsic change in 
the source. This is irrespective of whether the IDV results 
from intrinsic variability or from interstellar scintillation 
(ISS). If the variability is intrinsic, changes in the IDV 
timescale imply source size changes directly, while if the 
variability is due to ISS, it was argued that the change 
in timescale could be either the due to disappearance of 
the compact component or an increase in its angular size 



sufficient to quench the scintillation. Here we offer an al- 
ternative explanation based on an ISS origin of the vari- 
ability, showing that changes in the timescale need not 
reflect changes in source structure. 

Over the course of a year the motion of the Earth, 
projected on the plane of the sky is an ellipse with ec- 
centricity dependent on the ecliptic latitude of the source. 
Extragalactic sources show negligible proper motion, and 
hence the relative velocity of the ISM across the line of 
sight to the source is simply the difference between the 
ISM and the Earth's motion. If these are comparable in 
amplitude, then there may be a period when the vectors 
are parallel, and six months later, anti-parallel. In the first 
case the apparent ISM velocity as seen from the Earth is 
low, in the latter, high. 

For an ISS origin of IDV, it follows that at times 
of low relative velocity, the characteristic time-scale of 
the variability will be lengthened, while six months dis- 
placed, it will be correspondingly shortened. This form 
of behaviour has been seen in the rapidly variable IDV 
source J1819+3845 (Dennett-Thorpe & de Bruyn, 2000), 
and strongly supports the ISS origin for the IDV in this 
source. 

We have searched the literature on flux density mon- 
itoring of 0917+624 to determine its behaviour. The 
flrst observations were made in 1985 (Witzel et al. 1986; 
Heeschen et al. 1987) at 2.7 GHz with the Bonn 100 m 
telescope. Since then we have located a total of 10 well 
documented observing sessions undertaken with either the 
Bonn telescope or the VLA, and these are listed in Table 1. 

From this data we have determined the characteris- 
tic time scale, Tchar, (defined for convenience, as the mean 
time between successive peak-to-trough or trough-to-peak 
excursions) for each session at 2.7 and/or 5 GHz. This can 
only be done when the session was of sufficient length that 
several peaks/troughs are present in the data. Our choice 
of Tchar corresponds to the timescale on which the inten- 
sity structure function saturates, and hence closely follows 
the definition used by Kraus et al. (1999). Table 1 shows 
that Tchar is close to one day at 2.7 GHz, and roughly half 
this at 5 GHz. 
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Table 1 also shows that the observing sessions for 
0917+624 are not uniformly spread throughout the year, 
most are concentrated around the northern winter and 
spring. The only other session that we foimd that shows 
a slowing down of the IDV similar to that observed in 
September 1998 was the August 1985 session at 2.7 GHz 
(Heeschcn et al. 1987). These two observations 13 years 
apart clearly establish an annual cycle in this source. 

3. Effect of the Earth's orbital motion 

Changes in the variability timescale of a scintillating 
source occur due to changes in both the direction and 
speed of the scintillation pattern as it moves across the 
source-observer line of sight. The scintillation timescale is 



iiss = so{a)/v., 



app) 



(1) 



where SQ{a) is the length scale of the scintillation pat- 
tern in the direction parallel to the apparent scintillation 
velocity, a is the direction of motion relative to the pat- 
tern, and Wapp is the apparent speed of the scintillation 
pattern relative to the Earth. The scintillation timescale 
is formally defined as the timescale on which the auto- 
correlation function of intensity fluctuations reaches 1 /c 
of its saturation value; operationally, tiss ~ O.GTchar- We 
assume that the scattering occurs only at a single screen, 
so that the intrinsic velocity of the scattering medium rel- 
ative to the heliocentre is described by one velocity. 

3.1. Changes in the scintillation speed 

Annual changes in the scintillation speed result in an an- 
nual modulation of the variability timescale. The apparent 
scintillation speed is the magnitude of the component of 
the ISM's intrinsic velocity, relative to the Earth's veloc- 
ity, that is tangential to the line of sight. When the Earth's 
velocity nearly matches the intrinsic velocity of the scat- 
tering material, the scintillation pattern would then ap- 
pear to be stationary, causing the variability to stop. 

For a source at ecliptic latitude 6 and longitude (/), it 
is convenient to denote the Earth's velocity as 



sin ^f, cos "J", 0), 



(2) 



So (a) 



(4) 



where = 29.8 km/s is the Earth's orbital speed about 
the heliocentre and the orbital phase is 5* = 27rt — (j), and 
t is the time in years measured from the vernal equinox. 
The apparent scintillation speed is then 

t^app = [{v\\ +u®cos*)2 + K + t;©sin6'sin*)2]i/2, (3) Application to 0917+ 624. 



3.2. Changes in the scintillation direction 

Changes in the direction of the scintillation pattern, a, 
influence the scintillation timescale only if the scintilla- 
tion pattern is not circularly symmetric, i.e. if sq depends 
on a. Asymmetry in the scintillation pattern may be due 
to either intrinsic source structure or anisotropy in the 
scattering medium. 

The effect of source structure is important when the 
angular size of the source exceeds the angle to which a 
point source would be scatter-broadened by the scintilla- 
tion pattern. Opt- The scintillation lengthscale, sq, for a 
source of angular size 6 sip) is (e.g. Narayan 1992) 

Oe^t, 0sia)<9pt 
D9s{a), 9s{a)>ej,t' 

where D is the distance to the scattering screen. 

The angular scale, ^pt, depends on the strength of the 
scattering, which is a function of frequency. Above the 
transition frequency, i/q, the scattering is weak, and 9pt 
corresponds to the angular scale of the first Fresnel zone, 
Of = yc7(27fi7D). 

In the regime of strong scattering at frequencies below 
i/Q, one has 9pt = \J cj {2'kuoD){i' / vq)^'^ '^ . The exponent 
2.2 follows from the assumption that the turbulence in 
the scattering medium follows a Kolmogorov scaling (e.g. 
Armstrong, Rickett & Spangler 1995). The scale of size 
of the scintillation pattern of a point-like source increases 
strongly towards lower frequencies. 

IDV in extragalactic radio sources straddles the fre- 
quency range where weak and strong scattering occurs. 
This is particularly important in the case of 0917+624. 
Weak scattering is likely to be applicable to observations 
of 0917+624 at frequencies v ^ h GHz (Walker 1998), 
while the timescale of variability observed at 2.7 GHz is 
likely to correspond to refractive scintillation in the tran- 
sition to the strong scattering regime. 

Anisotropy in the scattering medium causes elongation 
of the scintillation pattern in a direction orthogonal to the 
magnetic field in the scattering medium (e.g. Goldreich 
& Sridhar 1995). Observations suggest elongations up to 
~ 2 : 1 are possible (e.g. Wilkinson, Narayan & Spencer 
1994; Spangler & Cordes 1998). The effect of anisotropy 
has been discussed by Backer & Chandran (submitted) 
on the scintillation parameters of nearby pulsars and IDV 
radio sources. Effects due to anisotropy become unimpor- 
tant for source sizes larger than O^i, since the scale of the 
scintillation pattern is then dominated by the angular size 
of the source itself. 



where U|| is defined as the speed of the scattering screen 
parallel to the Earth's ecliptic at closest approach, which 
occurs at ^ = 4>/2-k, and v±_ is its speed perpendicular 
to the ecliptic at this point, with positive values of v± 
indicating motion toward higher ecliptic latitudes. 



We show that the changes in the scintillation speed due 
to variations in the Earth's orbital velocity alone are 
sufficient to reproduce the behaviour of 0917+624. We 
consider only changes in the scintillation speed, and ne- 
glect any possible asymmetry in the source structure. 
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We also neglect the effect of anisotropy in the scatter- 
ing medium because the characteristic timescale of the 
intraday variability in 0917+624 is much longer than the 
Fresnel timescale, suggesting an angular size 0s (oi) ^ Op- 
(Numerically, one has 9p = 2/4 x 10^ m for a screen of 
£) = 200 pc at = 5 GHz.) 

The ecliptic latitude and longitude of this source are 
9 = 0.768 and (/) = 2.08 respectively. We find ISM ve- 
locities of = — 16 ± 3 km/s and = — 16 ± 3 km/s 
reproduce the sharpness of peaks observed Tchar at 2.7 
and 5 GHz. The scale size of the scintillation pattern, sq, 
may be calibrated by a measurement of the variability 
timescale given an estimate of the scintillation speed at 
that time. For the v± and v\\ required to reproduce the 
observed data this implies sq ~ 9.5 x 10^ m at 5 GHz and 
So « 1.9 X 10® m at 2.7 GHz. For a scattering screen at 
D = 200 pc, as used by Rickett et al. (1995), this cor- 
responds to angular sizes 6s = 30 and 65/uas at 5 and 
2.7 GHz respectively. This is quite consistent with the de- 
rived angular sizes of other IDV sources (e.g. Macquart et 
al. 2000). 

The comparison between the observational data and 
the model is shown in Figure 1. The model gives excel- 
lent agreement with the data at both 5 and 2.7 GHz. In 
particular, it clearly reproduces the peak in the August 
and September values of Tchar- The location of the peak 
observed in Fig. 1 is strongly sensitive to the ecliptic lon- 
gitude of the source, but only weakly dependent on the 
intrinsic velocity of the ISM. 

The validity of the above model can be readily checked 
with more data at both frequencies over the July through 
November period. Our model predicts a steady increase 
in the characteristic time scale during July, a significant 
slow-down in the variability through August and into 
September, with a decreasing time scale probably October 
through November. 

4. Discussion 

4-1- The intrinsic brightness temperature of the source 

Examination of the light curves during the slow-down pe- 
riod can also be used to set limits on any intrinsic vari- 
ability that may be present at radio wavelengths. During 
the slow-down the scintillations are virtually suppressed, 
so that any residual variations observed provide an upper 
limit to the intrinsic variations. This is an upper limit be- 
cause there are probably low-level scintillations remaining 
as the model of a simple screen at a single distance with 
a single velocity is certainly an oversimplification. 

For 0917+624, the August and September data still 
show the presence of low-level variability. Because of the 
~ 5 day extent of both sessions, the derived variability 
brightness temperature, estimated using a source diame- 
ter of 5 light days, is T^i, < 10^^ K, still well above the in- 
verse Compton limit, but significantly less than the 10^* K 



inferred if the intra-day variability were intrinsic. Fiitiirc 
observations of longer duration may better define the slow- 
down period and further reduce this limit. 

4-2. Earth Revolution Scintillation Synthesis 

The annual procession of the Earth around the sun and 
its effects on ISS shows a number of parallels with Earth 
rotation synthesis: in the latter we have a regularly con- 
structed interferometer that completes its synthesis in half 
a day (half of the Earth's rotation). In the former case, we 
have an irregular "interferometer" , namely the ISM, and it 
completes its "synthesis" in half a year, which corresponds 
to half of the Earth's rotation around the sun. 

We suggest that Earth Revolution Scintillation Syn- 
thesis can be used to model the microarcsecond struc- 
ture of compact radio sources. Once the ISM velocity has 
been determined for a particular source, the the shape of 
the scintillation light curves may provide information on 
source structure. 

The technique can be extended to studying the po- 
larization structure. Macquart et al. (2000) make use of 
this in studying the linearly and circularly polarized light 
curves for the strong and variable linear and circular polar- 
ization in the southern IDV source PKS 1519—273. Here 
both the linear and circular light curves follow the total 
intensity light curve quite closely, and it is clear that the 
microarcsecond source that is responsible for the IDV has 
the same structure in the total intensity, and linearly and 
circularly polarized components. 

5. Conclusion 

We have shown that a simple scintillation model accounts 
for the annual changes in the intraday variability timescale 
observed in the source 0917+624. These changes are due 
to the Earth's orbital motion, which influences the vari- 
ability timescale because its velocity with respect to the 
ISM changes with time. This causes a change in the speed 
at which the Earth moves across the scintillation pattern. 

This strongly suggests that ISS is the dominant mech- 
anism of intra-day variability in this source. 

We urge further observations of this source in the pe- 
riod July to November to test this model and constrain the 
velocity of the ISM transverse to the line of sight towards 
0917+624. 

We also suggest that observations of the source's vari- 
ability properties at various epochs can be used to gain in- 
formation on the two-dimensional structure in the source, 
as the direction of motion across the scintillation pattern 
varies during the year. 

Final Note During final preparation of this paper, 
we heard that Rickett, Witzel, Kraus, Krichbaum & 
Qian have independently discovered the annual cycle in 
0917+624 and the explanation in terms of ISS and are 
submitting an independent publication. 
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Figure 1: Comparison of the data at (a) 5 GHz and 
(b) 2 GHz with the scintillation model described in §3.3 
with W|| = —18 km/s, v± = —16 km/s. Scintillation pat- 
tern sizes of So — 9.5 x 10^ m and sq — 1.9 x 10^ m have 
been used at 5 and 2.7 GHz respectively (see text). 



